Relic keV sterile neutrinos and reionization 
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A sterile neutrino with mass of several keV can account for cosmological dark matter, as well as 
explain the observed velocities of pulsars. We show that X-rays produced by the decays of these relic 
sterile neutrinos can boost the production of molecular hydrogen, which can speed up the cooling 
of gas and the early star formation, which can, in turn, lead to a reionization of the universe at a 
high enough redshift to be consistent with the WMAP results. 
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There is an ample evidence that most of the matter in 
the universe is not ordinary atomic matter, but that it is 
made up of new, yet undiscovered particles. A seemingly 
unrelated astrophysical puzzle is the origin of the rapid 
motions of pulsars, whose velocities range in hundreds of 
km/s, while as many as 15% of pulsars have speeds in 
excess of 1000 km/s Q. It is intriguing that these two 
puzzles may have a simultaneous solution if there exists 
a sterile neutrino with mass of several keV and a small 
mixing with ordinary neutrinos. If such a particle exists, 
it would be produced in the early universe with the right 
abundance to be the dark matter 0, 0, , and it would 
have an important effect on the cosmological stucture 
formation |5j. The same particle would be emitted from 
a supernova with a sufficient anisotropy to give the pul- 
sar a kick velocity of the right mag nitude GIUEi. The 
pulsar kick from an anisotropic emission of sterile neu- 
trinos can enhance the convection during the first second 
of the supernova, which can increase the energy of the 
supernova shock and bring the supernova calculations in 
better agreement with observations Another hint at 
the existence of such a particle comes from the obser- 
vations of central galactic black holes with masses 
3 x 10 6 to 3 x 10 9 solar masses at redshifts as high as 
6.4. Dark matter in the form of sterile neutrinos could 
speed up the formation of supermassive black holes at 
high redshift 0] . Theoretical considerations lead one to 
believe that at least three sterile neutrinos should exist 
below the electroweak scale |ll|, in which case the neu- 
trino oscillations can explain the baryon asymmetry of 
the universe |l^.ll3|. 

Structure formation and reionization of the universe 
can be used to test some of the properties of dark-matter 
particles. Here we show that X-ray photons from the 
decays of the relic sterile neutrinos could help speed up 
the reionization of the intergalactic medium (IGM), in 
accordance with the WMAP observations of the early 
ionization at redshift z r = 17 ± 5 fl4l ll5T|. 

If the intergalactic medium (IGM) was reionized by 
stars (rather than some new physics 0), the WMAP 
data imply a very early and efficient star formation. Ster- 



ile neutrinos with keV masses have a non-negligible free- 
streaming length, which could delay the small-scale struc- 
ture formation 17], if it were not for the photons emit- 
ted in their decays. The width of the decay is extremely 
small, and the sterile neutrinos are stable on cosmolog- 
ical time scales. Nevertheless, some of them do decay 
producing photons with keV energies. If the dark matter 
is made up of sterile neutrinos, the X-ray photons from 
their decay can catalyze the formation of molecular hy- 
drogen in the amounts that can speed up the gas cooling 
and star formation, hence leading to an early reioniza- 
tion. 

The decays of particles during "dark ages" and their 
impact on reionization were discussed in Refs. 0, 0] . 
According to Mapelli and Ferrara > the electron frac- 
tion x e remains between 0.001 and 0.01 at redshifts 
10 < z < 100. These values are too small to produce the 
appreciable Thompson optical depth, but they are con- 
siderably higher than the electron fraction in the absence 
of sterile neutrinos. Mapelli and Ferrara concluded 
that the X-ray photons from sterile neutrino decays can 
play no role in reionization. However, even a small frac- 
tion of ions may be sufficient to catalyze the formation 
of H2, which, in turn, can precipitate a rapid star for- 
mation. Then the early stars can reionize the IGM. We 
will show that the electron fraction as high as 10~ 3 can, 
in fact, boost the production of molecular hydrogen, H 2 , 
to a value which is more than sufficient for a rapid star 
formation well before z x = 17 ± 5. 

Let us consider a singlet neutrino that has a non-zero 
mixing with the electron neutrino, and let us ignore other 
mixings for simplicity. Then the mass eigenstates have a 
simple expression in terms of the weak eigenstates: 

\v\) = cos 9 \v e ) — sin \v s ) (1) 
\v 2 ) = sin 6» |^ e ) +cos6»|i/ s ). (2) 

If the mixing angle is small, one of the mass eigen- 
states, v\ behaves very much like a pure v e , while the 
other, i>2, is practically "sterile" , which means it has weak 
interactions suppressed by a factor (sin 2 9) in the cross 
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section. The mixing enables the mass eigenstate v 2 to 
decay into lighter neutrinos, as well as v\ and a photon. 
The inverse width of v 2 — > i>i7 decay is 
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The X-ray photon produced in this two-body decay has 
energy E 7 — m s /2. 

The mass range consistent with v 2 being dark matter, 
with the pulsar kicks 0, and all the other constraints 

2 keV < m s < 8 keV. 

The corresponding photon energies are 

1 keV < E 1 < 4 keV. 

This range is rather conservative. A broader range of pa- 
rameters may be allowed, depending on cosmology |22( . 

The attenuation length A for real absorption of such 
photons in hydrogen [2ij] is 0.12 g/cm 2 < A < 7.7g/cm 2 . 
The dominant absorption process is the ionization edge 
of hydrogen [24| ■ Since the density of gas at redshift z is 
Po.h = 4 x 10~ 31 (1 + z) 3 gcm~ 3 , the energy attenuation 
length of a photons becomes smaller than the horizon 
size .ff -1 = Hq 1 ^ + z)~ z / 2 at redshift z > z t hick = 9 for 
E y = IkeV, and at z > z t hick = 148 for E 1 — 4keV. 
This includes a correction for the presence of helium. 
However, even at lower redshifts, z < zthick, the photons 
can ionize gas after they get redshifted. As the photon 
gets red-shifted, the cross section goes up with (1 + z)~ 3 
due to the ionization edge frequency behavior. The mean 
free path decreases as (1 + z)~ 3 / 2 . Hence, the photons 
emitted at redshift z a get all absorbed at some lower red- 
shift z b . For 4 keV photons, 7.7((1 + z b )/(l + z a )) 3 = 



0.0043(1 
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1 3/2 photons emitted at redshift 77 are 



all absorbed at redshift 40, implying a loss of about fac- 
tor 2 in the energy available for ionization and heating. 
Lower redshifts of the emission yield even greater losses 
in energy that can be used for ionization. One can safely 
assume that all the photons produced at redshift greater 
than z > 40 lose most of their energy to heating and ion- 
ization of the intergalactic medium within less than the 
Hubble time. 

Interactions of 1 — 4 keV photons with gas start out by 
ionization which absorbs the photon and ejects an ener- 
getic electron, which heats and ionizes further. Compton 
scattering transfers only a small fraction of the photon 
energy to the electrons. A cascade started by a single 
photon can ionize a number of atoms. Photoionization 
and Compton ionization of atoms are the dominant en- 
ergy loss processes [25| . The photon energy is divided be- 
tween reionization, heating, and excitations, with about 
30% of it going to ionization |2^|. Hence the number of 
ionized atoms can be estimated as ~ ry(m s /13.6eV) per 
X-ray photon, where r\ * 0.3. 
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FIG. 1: The rate k m of molecular hydrogen production as a 
function of redshift. For illustrative purposes, the gas tem- 
perature is assumed to be close to that of the CMBR (which 
is a good assumption, at least, for z > 150). 



The density of photons produced from the sterile neu- 
trino decays at time t is the sterile neutrino density times 
(t/r) — (to/r)(l + z)~ 3 / 2 . Each of these photons then 
ionizes ~ rj(ra s /13.6eV) atoms. In the absence of other 
effects, at redshift z>40, the fraction of ionized atoms is 
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Thus, the production of H + ions due to sterile neutrino 
decays occurs at a constant rate per atom, 



A = 1.4 x 10" 
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The free electrons are consumed in the reaction H + e — > 
#7 with the rate fei w 1.88 x 10- 10 (T K )- - 64 cm 3 s _1 , 
where Tk is the gas temperature in Kelvins. Hence, the 
electron fraction is described by the equation 



= A- kin H x z e = A -B I — 



1.6 



(0) 



where uh = ^o,ff(l + z ) 3 , n o,H = 2 x 10 7 cm 3 , and 
B = ki n ,H = 2 x 10- 17 s -1 . 

The star formation depends on cooling, which becomes 
very efficient when molecular hydrogen is produced. The 
critical Hi fraction needed for a rapid star formation is 
f c w 5 x lO" 4 [13. 

In the presence of ions, the molecular hydrogen can 
form in several reactions. One channel involves H ~ 



H 



H 



H~ +7 
H 2 + e- 



(7) 

(8) 



and proceeds at a high rate for z <C 200. However, since 
H~ is quickly destroyed at high temperature by the inter- 
actions with the cosmic microwave background (CMB) 
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FIG. 2: The fraction of ions, x e , in the absence of sterile 
neutrinos (thin line), and for the dark-matter sterile neutri- 
nos with masses 4 and 7 keV (lower and upper thick lines, 
respectively). Also shown is the limit of baryon decoupling 
from CMBR (dashed line). Below the dashed line, gas cooling 
is unaffected by the CMBR. 



photons, this channel is closed at higher redshifts. 
z 3> 200, the main reactions are 



+ H 
+ H 



Ht + 7, 
H 2 + H+ 



At 

(9) 
(10) 



Both channels depend on the presence of ions and free 
electrons, and, hence, on the value of x e . The fraction 
of molecular hydrogen can be found from the following 
equation 



/ = k m (t)n H (t)x e (t) (1 
« k m (t)n H {t)x e (t), 



x e (t)-2f(t)) 



(11) 



where k m is a function defined in Ref. 28] . The behavior 
of k m as a function of redshift is shown in Figure 1. It 
reaches a maximum fc„ l max rs 2 x 10 _16 cm 3 s _1 at z »s 
80. In the relevant range of parameters, the molecular 
hydrogen is not destroyed [27fl, and the right-hand side 
of eqn. (|11[) remains positive. 

Shortly after recombination, the ionization fraction 
x e r-j 1, and some molecular hydrogen is produced. How- 
ever, in the absence of sterile neutrinos, x e falls below 
10~ 3 by redshift z = 400, at which point further produc- 
tion of molecular hydrogen is stymied by the deficit of 
free electrons. In our case, the free electrons are contin- 
uously supplied by the decays of v s at the rate given by 
equation JSJ. 

We can now integrate equation © with the initial con- 
dition x e = 1 at recombination. The result can be used 
to obtain f(t) from equation 1|11|) . Figure shows the 
solution of equation 10 for three values of A correspond- 
ing to dark-matter sterile neutrino masses 4 and 7 keV, 
as well as for A = 0. 

Let us now estimate the fraction of molecular hydro- 
gen. Since H 2 is produced and not destroyed for z <C 10 3 , 



the value of f{z) is always greater than the fraction pro- 
duced at any earlier time. In particular, /(z) is greater 
than the fraction of hydrogen generated during one Hub- 
ble time at redshift z%, as long as z\ > z. In other words, 



f(z) > AJ( tl ) a J x t x 
k, k m (zi)n H {z 1 )x e {zi) 



(l + 2i) 3 / 2 ' 



(12) 



Let us consider a gas cloud that has virialized at some 
redshift z V ir- The evolution of this cloud takes place as 
usual ^2gj, as long as the electron fraction is not high 
enough to maintain the thermal equilibrium between the 
gas and the CMBR at redshift close to z v - lr . However, 
if the gas interacts strongly with the CMBR, the Comp- 
ton scattering keeps its temperature close to that of the 
CMBR, and the analyses of Ref. 28] do not apply. Comp- 
ton cooling is described by the following equation: 



dT 
It 
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where k, 
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= 2.6 x 10- 20 (1 + z) 4 



(13) 



To check the 



efficiency of the Compton cooling, one must compare the 
cooling rate with the expansion rate of the universe. This 
gives a limit on the fraction of ions. The Compton cooling 
is inefficient (and can be neglected) as long as 



x P < 0.9 x 10~ 3 
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This limit is shown in Figure as a dashed line. As 
one can see, for the parameters of interest, the decou- 
pling of gas from the CMBR occurs at z > 100. For 
smaller redshifts, one can use the semi-analytical results 
from Ref. 28] , making a correction for a higher electron 
fraction. 

One should also check that the gas temperature is not 
raised significantly by the presence of keV photons. If the 
gas temperature reaches 3000 K, the molecular hydrogen 
is destroyed. The total energy released in photons from 
the sterile neutrino decays is 



L = 



n Us (m s /2)r 1 
= 8 x 10- 30 ^- 
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One can compare this heating process with the cooling 
processes discussed in Ref. |28| ■ For the relevant range of 
parameters, the heat deposited by sterile neutrinos can 
be neglected. 

Let us now consider a gas cloud collapsing and viri- 
alizing at some redshift z v j r . The density in a recently 
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virialized cloud is about a factor 187r 2 higher than the 
mean hydrogen density |28| : 



23 cm" 3 



l+£ 
100 



Let us consider z v ; r = 100. The fraction of molecular hy- 
drogen produced at this redshift, according to eqn. l|T2*)l. 



/ > Aioo/ = 9 x 10- 4 

km. \ f X \ / 

71 J VlF 3 / V 



20 cm- 



(16) 



This lower bound on the molecular hydrogen fraction 
is higher than the critical value f c « 5 x 10~ 4 derived 
by Tegmark et al. for a successful collapse, especially 
near the maximum of k m . This maximum is achieved at 
z — 80, as shown in Figure Q an d one expects the most 
efficient production of Hi for z v - lr ~ 10 2 . In the absence of 
sterile neutrinos, the electron fraction is around 3 x 10~ 4 , 
for which / drops below f c . 

In conclusion, if sterile neutrinos with masses of sev- 
eral keV constitute the dark matter, their decays pro- 
duce enough photons to boost the production of molec- 
ular hydrogen in clouds collapsing at redshifts as high 
as z ~ 100. The presence of molecular hydrogen facili- 
tates rapid cooling and the early star formation. At the 
same time, the sterile neutrinos can have a non-negligible 
free-streaming length, which can affect the early halo col- 
lapse 0, El • The full analysis must combine the enhanced 
production of molecular hydrogen with simulations of 
dark matter halo formation. Such detailed analysis is 
beyond the scope of this Letter; it will be presented else- 
where. 
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